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LONG-TERM  GOALS 

The  long-term  goals  of  this  project  are  to  quantify  turbulence  and  to  understand  the  dynamics  and 
implications  of  turbulent  mixing  in  the  coastal  ocean. 

OBJECTIVES 

The  specific  objectives  of  this  project  are  to  obtain  a  critical  evaluation  of  turbulence  closure  models 
that  are  commonly  used  in  the  coastal  ocean,  and  to  identify  ways  to  improve  these  models.  Of 
particular  interest  are  the  MY-2.5  model  (Mellor  and  Yamada,  1974,  1982)  and  the  k-epsilon  model 
(e.g.,  Rodi,  1987),  which  involve  coupled  equations  describing  turbulent  kinetic  energy  and  turbulent 
scale.  Also  of  interest  is  the  KPP  model  (Large  et  al.  1994),  which  determines  the  eddy  viscosity  and 
diffusivity  in  terms  of  boundary  layer  thickness  and  boundary  fluxes  of  momentum  and  buoyancy.  In 
the  present  applications,  work  to  date  has  clarified  the  turbulence  energetics  (see  below),  so  that  the 
primary  focus  is  on  processes  controlling  turbulence  scale,  which  is  by  far  the  most  uncertain 
component  of  turbulence  closure  models. 

APPROACH 

The  approach  is  to  compare  computations  based  on  turbulence  closure  models  with  oceanic  turbulence 
measurements  obtained  in  ONR's  Coastal  Mixing  and  Optics  (CMO)  program  and  in  an  NSF-funded 
study  in  the  Hudson  estuary.  Turbulence  closure  models  produce  estimates  of  turbulent  fluxes,  given 
the  Reynolds-averaged  profiles  of  velocity  and  density.  In  the  CMO  and  Hudson  observational 
studies,  both  turbulent  fluxes  and  Reynolds-averaged  profiles  were  measured.  Thus  the  measurements 
permit  direct  tests  of  turbulence  models,  in  contrast  to  the  more  common  indirect  tests,  based  on 
measurements  of  Reynolds-averaged  profiles  and  the  assumption  of  one-dimensional  dynamics  to  infer 
turbulent  fluxes.  In  addition  to  measurements  of  turbulent  fluxes,  the  CMO  and  Hudson  data  sets 
include  measurements  of  dissipation  rates,  which  permits  diagnosis  of  the  turbulence  energetics. 

WORK  COMPLETED 

Work  on  this  project  has  focused  on  two  tasks.  The  first  is  comparison  of  computations  based  on 
turbulence  closure  models  with  measurements  obtained  in  the  CMO  program.  Some  of  this  work  is 
described  by  Shaw  and  Trowbridge  (submitted).  The  second  task  is  analysis  of  a  unique  set  of 
measurements  obtained  from  a  spatial  array  of  near-bottom  velocity  sensors,  which  were  deployed  in 
shallow  water  near  Scripps  pier,  in  California.  Although  not  described  in  the  original  proposal  for  this 
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project,  the  California  measurements  were  pursued  because  they  explicitly  determine  turbulence  scale, 
which  must  be  inferred  from  temporal  variability  in  conventional  turbulence  measurements.  The 
analysis  of  the  California  measurements  is  described  by  Trowbridge  and  Elgar  (submitted). 

RESULTS 

Results  of  the  analysis  of  CMO  measurements  include  successful  closure  of  simplified  balances  for 
turbulent  kinetic  energy  and  momentum  (Figure  1)  and  an  evaluation  of  MO  similarity  theory  (a 
simple  turbulence  closure  model)  and  the  KPP  model  (Figure  2).  Detailed  explanations  of  these  results 
are  provided  in  the  figure  captions.  Of  particular  interest  is  failure  of  MO  similarity  theory  (Figure 
2a),  in  contrast  to  crude  but  quantative  agreement  between  measurements  and  computations  based  on 
the  KPP  model  (Figure  2b),  because  of  incorporation  of  the  effect  of  finite  boundary  layer  thickness  on 
turbulence  scale. 

Results  of  the  analysis  of  the  California  measurements  includes  explicit  determination  of  the  scales 
contributing  to  the  dynamically  important  alongshore  component  of  the  turbulent  Reynolds  shear  stress 
(Figure  2).  Of  particular  interest  is  demonstration  that  turbulence  scales  shorten  during  upwelling- 
favorable  alongshore  flows,  probably  because  of  the  intensification  of  near-bottom  stratification,  which 
limits  turbulence  scales. 

IMPACT/APPLICATIONS 

This  work  will  produce  critical  tests  and  ultimately  improvements  of  turbulence  closure  models  that 
are  widely  used  in  the  coastal  ocean. 

TRANSITIONS 

None. 

RELATED  PROJECTS 

Trowbridge's  and  Y.  C.  Agrawal's  participation  in  the  ONR  HYCODE  program  are  capitalizing  on  the 
insights  about  turbulence  dynamics  and  measurements  that  are  being  obtained  as  part  of  this  study. 
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Figure  1:  Tests  of  (a)  a  turbulent  kinetic  energy  balance  and  (b)  an  alongshore  momentum  balance, 
based  on  measurements  over  the  New  England  shelf  during  1996  and  1997  as  part  of  ONR's  Coastal 
Mixing  and  Optics  (CMO)  program.  The  measurements  were  obtained  with  an  array  of  BASS 
acoustic  travel-time  current  meters  ( Williams  et  al.,  1987)  at  heights  between  0.4  and  7.0  m  above 
bottom.  The  energy  balance  tested  in  (a)  is  shear  production  (P)  equals  dissipation  ( (c)  plus 
buoyancy  flux  (B).  In  this  example,  shear  production  and  dissipation  are  much  larger  than 
buoyancy  flux.  The  balance  P  =  B  +  s  closes  well.  The  momentum  balance  tested  in  (b)  is 
integrated  over  the  thickness  of  the  bottom  boundary  layer,  and  involves  temporal  acceleration, 
Coriolis  acceleration,  pressure  gradient,  and  bottom  stress.  In  (b),  p0  is  a  fixed  reference  density,  f 
is  the  Coriolis  parameter,  z  is  height  above  bottom,  t  is  time,  v  is  cross-shelf  velocity,  subscript 
8  denotes  evaluation  at  z  =  8 ,  where  8  is  boundary  layer  thickness,  and  rbx  is  the  alongshore 
component  of  the  turbulent  Reynolds  shear  stress,  extrapolated  to  the  sea  floor.  This  calculation  is 
based  on  8-40  m.  The  momentum  balance  closes  well.  In  both  panels,  r 2  is  the  squared 
correlation  coefficient  and  b  is  the  regression  coefficient  for  the  left  and  right  sides  of  the  equation. 
Successful  closure  of  the  turbulent  energy  and  momentum  balances  implies  that  the  measurements 
of  stress,  dissipation,  and  buoyancy  flux  are  accurate  and  that  the  turbulence  energetics  and  bottom 

boundary  layer  dynamics  are  relatively  simple. 
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Tot  of  MO  model  at  z  =  1 ,3  m  (r2  =  0  64.  b  =  0.22±0.0t ) 


Test  of  KPP  model  at  z  -  1.3  m>  (rz  -  0.6G.  b  ■  0.56+0.03) 


Figure  2:  Tests  of  turbulence  closure  parameterizations,  including  (a)  Monin-Obukhov  (MO) 
similarity  theory  and  (b)  the  K  Profile  Parameterization  proposed  by  Large  et  al.  (1994),  based  on 
measurements  obtained  during  the  CMO  experiment,  z  is  height  above  bottom,  u  is  the  Reynolds- 
averaged  horizontal  velocity  vector,  r,t  is  the  alongshore  component  of  the  turbulent  Reynolds  shear 

stress,  8  is  boundary  layer  thickness,  r2  is  squared  correlation  coefficient,  and  b  is  regression 
coefficient.  In  both  panels,  the  Reynolds-averaged  velocity  and  density  profiles  are  inputs  to  the 
model,  zzx  is  the  model  output,  and  model  estimates  of  zzx  are  compared  with  estimates  obtained 

directly  from  the  turbulence  measurements.  In  the  MO  model,  zzx  depends  only  on  p0,z,d/  Gz,  and 
9dp  /  oz ,  where  p  is  density  and  8  is  gravitational  acceleration.  In  the  KPP  model,  z7X  depends,  in 

addition,  on  the  thickness  of  the  boundary  layer.  Both  MO  and  KPP  estimates  of  stress  are  well 
correlated  with  measured  stresses.  The  MO  estimate  is  smaller  by  a  factor  of  five  than  the 
measurements.  The  KPP  estimate  is  better  because  it  incorporating  dependence  on  boundary  layer 

thickness. 
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Shape  of  structure  lurvction  tor  downweJIing- favorable  alongshore  flow 


Figure  3:  Estimates  of  the  quadratic  structure  function  Dvc0  (j  y),  obtained  from  measurements  at 

a  water  depth  of  approximately  4  m  near  Scripps  pier,  in  California.  In  this  application,  the 
structure  function  describes  the  scales  contributing  to  the  alongshore  component  of  the  turbulent 
Reynolds  shear  stress.  The  measurements  were  obtained  from  an  alongshore  array  of  velocity 
sensors.  0  y  is  alongshore  separation,  and  D  jrnax  is  the  maximum  separation  resolved  by  the  array. 

The  measurements  are  compared  with  an  expression  based  on  measurements  in  the  neutrally 
stratified  atmospheric  surface  layer,  in  which  the  ratio  of  D  y)  to  Dua)(j~  >max)  depends  only  on 
1  v  /  z ,  where  z  is  height  above  bottom.  During  downwelling-favorable  alongshore  flows,  the 
measured  structure  function  is  consistent  with  the  atmospheric  expression.  During  upwelling- 
favorable  alongshore  flows,  however,  turbulence  scales  are  shorter  than  indicated  by  the 
atmospheric  expression,  probably  because  of  near-bottom,  stable  stratification,  which  limits 
turbulence  scale  and  intensifies  during  upwelling.  Turbulence  measurements  from  spatial  array  of 
sensors  are  nearly  unique  in  the  ocean,  as  is  the  explicit  demonstration  that  turbulence  scales 

shorten  during  upwelling. 
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